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Abstract 

We examine the effects of boundary kinetic terms in the Randall-Sundrum model 
with gauge fields in the bulk. We derive the resulting gauge Kaluza-Klein (KK) state 
wavefunctions and their corresponding masses, as well as the KK gauge field couplings 
to boundary fermions, and find that they are modified in the presence of the bound- 
ary terms. In particular, for natural choices of the parameters, these fermionic cou- 
plings can be substantially suppressed compared to those in the conventional Randall- 
Sundrum scenario. This results in a significant relaxation of the bound on the lightest 
gauge KK mass obtained from precision electroweak data; we demonstrate that this 
bound can be as low as a few hundred GeV. Due to the relationship between the 
lightest gauge KK state and the electroweak scale in this model, this weakened con- 
straint allows for the electroweak scale to be near a TeV in this minimal extension of 
the Randall-Sundrum model with bulk gauge fields, as opposed to the conventional 
scenario. 
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1 Introduction 



The Randall- Sundrum (RS) model pQ offers a new approach to the hierarchy problem. 
Within this scenario, the disparity between the electroweak and Planck scales is generated 
by the curvature of a 5-dimensional (5-d) background geometry which is a slice of anti-de 
Sitter (AdSs) spacetime. This slice is bounded by two 3-branes of equal and opposite ten- 
sions sitting at the fixed points of an S l /Z2 orbifold which are located at <fi = 0, ir where <fi is 
the coordinate of the 5 th dimension. The 5-d warped geometry then induces an effective 4-d 
scale A w of order a TeV on the brane at <fi = it, the so-called TeV brane. A„- is exponentially 
smaller than the effective scale of gravity given by the reduced Planck scale, Mpi, with the 
suppression being determined by the product of the 5-d curvature parameter k and the sep- 
aration r c of the two branes, i.e., = Mpie~ kVcn . In this theory, the original parameters 
of the 5-d action are all naturally of the size ~ Mpi, while in the 4-d picture a hierarchy 
with A n /Mpi < 10 -15 appears. It has been demonstrated [2j that this scenario is naturally 
stabilized without the introduction of fine-tuning for kr c ~ 11, which is the numerical value 
required to generate the hierarchy. This model leads to a distinct set of phenomenological 
signatures that may soon be revealed by experiments at the TeV scale; these have been 
examined in much detail 0IUI2I- 

For model building purposes, it is advantageous to extend the original framework of 
the RS model, where gravity alone propagates in the extra dimension, to the case where at 
least some of the Standard Model (SM) fields are present in the bulk [6 . The simplest such 
possibility is to have the SM gauge fields in the bulk jHE], while the other SM fields remain 
localized on the TeV brane. This scenario was examined in detail some time ago |1] where it 
was found that the couplings of the Kaluza-Klein (KK) excitations of the bulk gauge fields 
to matter on the boundary were enhanced by a factor of y/2nkr c compared to those of the 
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corresponding zero- mode states. Since kr c ~ 11 as discussed above, this enhancement is 
numerically significant and is approximately a factor of ~ 8.4. Severe constraints are then 
imposed on the gauge KK excitations from precision electroweak measurements; these imply 
that the mass of the lightest gauge KK state must be heavy with mi > 25 — 30 TeV. Due to 
the relation in this model between the masses of the KK states and the scale of physics on 
the TeV brane, this bound then correspondingly requires that > 100 TeV. This places 
such a scenario in a very unfavorable light in terms of its original motivation of resolving the 
hierarchy. 

Here, we investigate whether Brane Localized Kinetic Terms (BLKT's) for bulk gauge 
fields can modify these results. Recently, Carena, Tait and Wagner |7| examined the phe- 
nomenological influence of such localized kinetic terms for bulk gauge fields within the con- 
text of flat TeV _1 -sized extra dimensions. These authors showed that such terms can lead to 
significant modifications to the gauge KK spectrum, as well as to the KK state self-couplings 
and couplings to fields remaining on orbifold fixed points. Localized kinetic terms are ex- 
pected to be present on rather general grounds in any orbifolded theory, as was shown by 
Georgi, Grant and Hailu jH], since they are needed to provide counterterms for divergences 
that are generated at one loop. In addition, these boundary terms can have important im- 
plications in a number of different situations such as model building jU], the construction 
of GUTs in higher dimensions with flat ^U] or warped ^T] geometries, and, in a different 
context, the quasilocalization of gravity [T2]. Up to now the phenomenological implications 
of BLKT's in the case of warped extra dimensions have not been examined. 

We examine the effects of localized kinetic terms in the case where the SM gauge fields 
are present in the RS bulk and the remaining SM fields are on the TeV-brane. As in the 
case of TeV _1 -sized extra dimensions, we find that these boundary terms alter both the KK 
spectrum as well as their associated couplings to the remaining wall fields, but in a manner 
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which is quite distinctive due to the curved geometry in this scenario. We will show that 
for a reasonable range of model parameters these modifications can be quite significant and 
lead to a natural reduction in the strength of the KK couplings. This allows the restrictions 
on A n to be reduced by more than an order of magnitude and it would then no longer be 
unfavorable to have gauge fields be the only SM fields present in the RS bulk. The next 
section contains a detailed discussion of the modification of the usual gauge boson analysis 
within the RS framework due to the existence of brane localized kinetic terms. In Section 
3 we discuss the implication of these results for the phenomenology of this scenario and our 
conclusions can be found in Section 4. Special limits of the KK mass equation are discussed 
in an Appendix. 

2 Formalism 

In this section, we derive the wavefunctions and the KK mass eigenvalue equation for bulk 
gauge fields with BLKT's in the RS scenario. We then compute the coupling of the gauge 
field KK modes to 4-d fermions localized on the TeV brane, where the scale of physics is 
of order the electroweak scale. Our derivation will demonstrate how the boundary terms 
modify these couplings from their original form and result in a possible softening of the 
precision electroweak bounds on the mass of the lightest KK gauge state. In the rest of the 
paper, KK modes refer to those from the bulk gauge fields, unless otherwise specified. 

We perform our calculations for a U(l) gauge field; the generalization to the case of 
non-Abelian fields is straightforward. Let us start with the 5-d action for the gauge sector 




+ 




(1) 
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where the terms proportional to cq and represent the kinetic terms for the branes located 
at = 0, it, respectively, and the metric Gmn is given by the RS line element 

ds 2 = e- 2u T]^dx^dx u - r 2 c d<p 2 ; a = kr c \<p\ . (2) 

Here the dimensionless constants Co j7r would naturally be expected to be of order unity. Note 
that we can neglect terms proportional to d^A^ on either brane as demonstrated by [7]. We 
will work in the A§ = gauge so that terms including A§ can also be neglected. Thus only 
the usual 4-d kinetic term appears in the brane contributions to the action. In our notation, 
A,B,... = 0, . . . , 3, 5 and a, ft, . . . = 0, . . . , 3. Here, k is the 5-d curvature scale, r c is the 
compactification radius, and x 5 = r c (p with — n < (f> < ir. The field strength is given by 

Fmn = 9m An — djs[A M . (3) 
The gauge field is assumed to have the KK expansion 

^(^) = EfW^' ( 4 ) 

n V ' c 

Substituting the above expansion in the action Q yields 

Sa= ~ ^~ fd'xY, {[l + c 5(<P) + c n 5(<P-rr)]F^F^x {m) X in) 

^ r c J m,n 

- (2/r c 2 )^ m )4") e - 2CT [( C //rf0)x (m) ][(rf/rf0)x (n) ]} , (5) 

where now all the Lorentz indices are contracted using the 4-d Minkowski metric. In order 
to cast this action in the diagonal form, we demand 

j dcf> [1 + c 5(0) + c n 5{<P - tt)] x (m) X {n) = ZJ mn (6) 
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and 
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(m) 
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2 xmn 



These equations imply 



^ ( 6 " 2CT ^ x(n) ) + [1 + Co ^ + ^ 5(0 ~~ n)] r ° m " X(n) = ° 



Away from the boundaries at = 0, ir, the solution to Eq.(JE} is given by jU IHj 



(7) 



X W (0) = T r Cl(^) 



(9) 



where = {m n /k)e a ; ( q (z n ) = \J q (z n ) + a n Y q (z n )), with J q and F g denoting the usual 

Bessel functions of order q. The normalization N n is set by the condition 
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implying that the zero mode is given by x^ = l/\/27r. Integrating Eq.flSJ) around the fixed 
points = and <fi = ±7r yields 



d 



2—X {n \0)+c r 2 c m 2 n x {n) (0) = 



:ir. 



and 



24 (n) (^)-e 2 ^c,r c 2 m ^ (n) (^)=0 
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respectively. 

Let e n = 2 n (0). For light KK modes with masses m n ~ 1 TeV, which are of phe- 
nomenological interest, e n <C 1; for k ~ Mp; we have e n ~ 10~ 15 . Substituting for x^ from 



Eq.® in Eq.lfTTj). and using the identity 



(13) 



we obtain 



Jo( £ n) + &0£nJl{£n) 



(14) 



io(e«) + fe^l^n) 



where #o = cokr c /2 and a prime (') denotes a derivative. Note that with Co of order unity, 
we may expect values of 5q of order 10 since kr c takes on the value ~ 11 in order to address 
the hierarchy. Similarly, from Eq. (jl2j) . we obtain 



where 5 n = c w kr c /2 and x n = z n (ir). The roots x n of Eq. ffTo]) yield the gauge KK mass spec- 
trum, m n = x n ke~ kTc7T . We will study various limits of this root equation in the Appendix. 
Following the same arguments as above, we may expect that |5 07r | < O(10), assuming no 
fine-tuning. 

Due to the exponentially small value of e n , it may first appear that the (^-dependent 
corrections to a n are highly suppressed, and that a n is small, of order ~ 0.01, as in the 
standard RS scenario. However, an expansion of Eq. (pPf]) in the parameter e n yields 



where 7 ~ 0.577 is Euler's constant. a n remains small over most regions of parameter space 
of interest here, except for when 5q ~ —kr c i\ ~ —35.4. There are two sources of modifications 
to a n as compared to the original RS framework due to the presence of the boundary terms: 
(1) the appearance of the term proportional to 5q in the denominator which is not suppressed 
by the exponential warp factor, and (2) the BLKT-modified values of x n . In our numerical 



( {x n ) - SnXndtXn) = 



(15) 



vr/2 



(16) 



ln(x„/2) - kr c ir + 7 - 5 ' 
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Figure 1: (a) The behavior of the first five roots as 5 n is varied, assuming So = 0. Note 
the sharp increase in x n near 5 n = for all KK levels, (b) The large negative 5^ region is 
expanded to show the new root originating at the value 5 n = —kr c n ~ —35.4. 
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analyses we will use the exact expression for a n , rather than the approximate version given 
above, and we take kr c = 11.27. 

To get a flavor of how the KK tower mass spectrum is modified in this scenario with 
BLKT's, we must solve Eq. lfTKj) for a range of values of 5o j7r . We first take the simplest 
possibility by setting 5q = 0, and examine the roots function of 5^. Our results are 

displayed in Fig. ^ From Fig. we see that in the region near 5 n = the roots undergo a 
substantial shift but remain relatively independent of variations in S„ away from the origin. 
For large negative values of S n , the roots are almost the same as those obtained in the case 
of the graviton KK spectrum in the original RS framework |3] . This can be seen analytically 
by taking the limit of large \5 n \ in Eq.lfTHjl and neglecting a„Cl. When the value of — 5 n 
exceeds kr c ir, a new root appears which is barely perceptible in Fig. [TJi, but whose turn-on 
is shown in detail in Fig. ^p. The analytical source of this root is discussed in the Appendix. 
We will later see that this signals the appearance of an unphysical ghost, i.e., a tachyonic 
state with negative norm. 

When 5q ^ the modifications to the roots are sensitive to the sign of 8q. For 
$o — $it, and in the range 5 n > —30, the roots appear to be very similar to those obtained 
when 5q = 0, as is shown in Fig. |2^l. The major difference in this case is the location of 
where the additional ghost root materializes; as shown in the Appendix, this root appears 
at the value 5 W = —(kr c n + 8 ). When 5 > (<)0 the turn-on of this root thus moves to 
smaller (larger) values of S n . In the case 8q = S n , Fig. |2k shows that the new root appears 
at 5 n = —kr c 7r/2 as expected. We also see that for 5 n ^ — 7ckr c the roots approach those 
which govern the KK graviton spectrum, as discussed above. In Fig. |2b, we display the 
roots for the case S = —5 n . For most values of S n shown in the figure, these are found to be 
similar to the case where 5q = except that a new light root does not appear in this case. 
This is as expected, based on the discussion in the Appendix. When 5 n is large and positive, 
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<5tt ^ vrA;r c , the roots again approach those for the case of KK gravitons, as anticipated. 

When 5 n = 0, the roots exhibit a weak dependence on 5 away from the 5q = —kr c 7r 
region; this can be seen in Fig. EJ The 5q dependence of the roots is isolated in the coefficients 
a n as seen in Eq.()16j). From this we see that the roots have a slightly negative slope with 
increasing S . 

3 Analysis 

In this section, we study the couplings of the KK modes to the boundary fermions at <fi = tt. 
The values of these couplings are the driving force behind the important constraints on 
the mass of the lightest KK gauge field in the original RS scenario [I]. In this section, we 
reexamine these bounds in the presence of brane localized kinetic terms. To begin, we note 
that the diagonalization conditions (JUJ) and (J7J) yield the 4-d action 

Sa = -\I d 4 x^2 Z n (F^F$ - 2m 2 n A^A^) , (17) 

where we have 

Z n = l + c [x (n) (0)] 2 + c.[x (n) (vr)] 2 . (18) 

Note that for physical fields, we demand that Z n > for all n > 0. 

The coupling of a bulk gauge field to fermions ip localized at (p = it is given by 
the action 

S^ = - J d(j)d*x{detV)g h $VZ 1 a A ll {x 1 <j ) )^8( ( l)->K)i (19) 
where V£ = e a rj£ is the vielbein, detV = e~ 4cr , and g 5 is the 5-d coupling constant. Using 
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Figure 3: The behavior of the first five roots as a function of 5q taking 5^ = 0. 



the expansion (J2J) in Eq. (fTT)|) . we obtain 



S,i, = — d x 



95 



„7rfcr t 



(20) 



where we have performed the redefinition ip — > e 3fcr W 2 ^ t m ake the ip kinetic terms canon- 



ical. Here, we have used the normalization 



Nt 



x 



n /-2 



kr c e 2 n 



(t(x n )W. 



(21) 



where 



W = (1 - 28* + 5y n ) 



1 - 



e 2 n e i (e n )[l + 25 + el5" ] 
xlCl{x n )[l-25 v + xl5l} 



(22) 



Note that this normalization is mode dependent. To bring the gauge field kinetic terms in 
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Eq. (jl7|) into the canonical form, we require 



(23) 



Eqs.flHJ) and (J21 then yield 



S,i, — — d x 



g jrfA®il> + E ^7 M 4 n V 



(24) 



with 



6-5 



y/2nr c Z 



(25) 



and 



5- n = g J2nkr i 



WZ„ 



(26) 



Here, is identified as the usual 4-d coupling constant for the interactions of the fermions 
with the zero- mode KK state. The KK excitation couplings are given by g n , which due 
to the presence of the boundary terms are now mode dependent. Without the presence 
of the brane kinetic terms, one obtains Z = Z n = 1 which yields the result that g n is 



an approximately mode-independent constant given by g n ~ go\ / 2irkr c , since in this case 
a n ~ 10~ 2 and W ~ 1. 

We can now compute the size of the effect of the brane localized kinetic terms on the 
ratio of couplings g n / go- Using Eq. (fTH|) . we obtain 



1 + 



£p + £tt 

7rkr r 



(27) 



and, neglecting terms of order e n , 



Z ~ 



1 + W' 1 



25-n 



(28) 
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Figure 4: The ratios of the first five KK couplings to that of the zero-mode as a function of 
5-x for the cases (a) 5 = 0, (b) 5o = 5 n , and (c) 5o = —S n . The lightest (heaviest) KK mode 
corresponds to the top (bottom) curve in each case. 
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Note that similar to the case of TeV _1 -sized extra dimensions, Z n is found to be n-dependent. 
From the requirement that Z > 0, we obtain the constraint Sq + 5 n > —irkr c ~ —35.4. If 
this condition isn't satisfied, the state then has negative norm and is unphysical. 

We first consider the case S n > 0, in order to avoid the ghost states. Figures HK-c 
display the existence of a significant fall-off of the ratio of couplings g n /go as a function of 
S n . The values of the coupling ratio in the negative S„ region, away from the regions where 
the ghost states appear, are shown in Fig. |^l-c. In both cases we see the mode dependence 
of the coupling strength away from the value S n = 0. Note that the coupling strength of the 
first KK excitation is substantially smaller when S n < compared to the case with positive 
values of S n . It is clear that compared to the original RS framework, the ratio g n / g is now 
naturally much reduced. 

Due to the shifts in the mass spectrum and the reduced values of the ratio g n j g , 
we now expect that the inclusion of the localized kinetic terms can result in a significant 
loosening of the precision electroweak bounds jlj on the RS model. To obtain the constraints 
on the KK gauge boson mass spectrum, we need to consider the quantity 

V iSl ,S a )^{ S -f{-) 2 , (29) 

where the x n are the roots given above, and examine its detailed dependence on S 0>n . 
This quantity describes the relative strength of the set of contact interactions induced by 
dimension-6 operators arising from KK gauge exchanges in precision electroweak observ- 
ables [4J; a full description of the analysis procedure and the dependence of the precision 
electroweak variables on V can be found in Ref. ^Hl- Here, we perform a global fit to 
the most recent precision electroweak data sample with the use of ZFITTER. Since 
for the case Sq^ = we know from our previous analysis that precision measurements 
imply the bound mi\ns ^ 25 TeV for the lightest gauge KK excitation jlj, we can eas- 
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Figure 5: The ratios of the first five KK couplings to that of the zero-mode as a function of 
negative 5 n for the cases (a) So = 0, (b) So = S n , and (c) So = —S n . The lightest (heaviest) 
KK mode corresponds to the top (bottom) curve in each case. 
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ily obtain the corresponding constraint for non-zero values of <5o i7r . We obtain the bound 

ttiiIblkt > rn\\Rs\y 8q)/V(0, 0)] 1 / 2 ; the numerical results of this analysis is presented in 
Fig. IHk-b. Here, we see that the bound on m\ falls rapidly as the magnitude of 8^ increases. 
For example, with 8 n ~ 10 — 20, we see that mi can be as low as 2 — 3 TeV which will 
lead to an observable signal in Drell-Yan production for the first gauge KK excitation at the 
LHC This is demonstrated in Fig[7|for the case of mi = 5 TeV. For larger values of 8 n , such 
states may be observable at the Tevatron. For 8 n < away from the ghost root region, the 
constraints from electroweak data are even weaker; here, direct searches at the Tevatron can 
place bounds on the parameter space. 

So far, we have shown that the inclusion of boundary kinetic terms allow for a sub- 
stantially lighter mass spectrum for the gauge KK states than in the conventional RS sce- 
nario. However, the question remains as to whether this ameliorates the stringent bound of 
A,r > 100 GeV. To address this issue, we convert the constraint on mi into one on A„- by 

noting that m 1 = xiA w k/ Mpi, recalling that k/Mpi < 0.1, and using the values of x\ that 
we computed above. Solving for A,,-, we find the results displayed in Fig. |H1 Note that for 
most values of the parameters, the resulting bound on A„- remains quite large for positive 
values of 8 n . This is due to the interplay of the modified roots and the constraints on the 
gauge KK mass spectrum. However, for 8 n < 0, we see that the lower bound on A,,- is easily 
reduced by almost an order of magnitude or more. This would seem to be the preferred 
region of parameter space for this model in its relation to the hierarchy problem. 

4 Conclusions 

A minimal, yet interesting, extension of the RS model is obtained by assuming that only 
gauge fields propagate in the bulk. However, this setup has been shown to be constrained by 
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Figure 6: (a) Lower bound on the mass of the lightest gauge boson KK excitation as a func- 
tion of 8-x > from a fit to the precision electroweak data for the cases 8 = 5 T (top/green), 
5 = (middle/red), and 5 = —8^ (bottom/blue), (b) Same as above but now for the region 
S w < with 5 = 5 W (bottom/green), <5 = (middle/red), and S = —8^ (top/blue). 
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M (GeV) 

Figure 7: Drell-Yan event rate at the LHC for the first gauge KK state in the RS scenario 
with boundary terms. We have assumed a luminosity of 300 fb _1 and have set mi = 5 TeV. 
The steeply falling histogram represents the SM Drell-Yan continuum. 

precision electroweak data to have a scale A n ^ 100 TeV for physics on the TeV brane jl]. 
This makes such a scenario less attractive as a means for resolving the gauge hierarchy. The 
stringent bound on the electroweak scale A n is a result of the strong coupling of the gauge 
KK modes to the boundary fermions. It is thus interesting to examine if a modification of 
the boundary physics can lead to a relaxation of this constraint. One such possibility which 
is a simple extension of the theory, motivated by field theoretic considerations (SJ, is the 
addition of brane localized kinetic terms for the bulk gauge fields. 

In this paper, we assumed the presence of localized gauge field kinetic terms on both 
the Planck and TeV branes in the RS model. We then derived the wavefunctions and the 
mass spectrum of the gauge KK modes. The couplings of the KK states to the SM brane 

localized fermions were then calculated. It was shown that for natural choices of parameters, 
a substantial suppression of these KK couplings, compared to the original model, can be 
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Figure 8: Bounds on A n as a function of for the cases 5q = 5 n (green), 5q = (red), and 
5q = —5^ (blue), as derived from the constraints on the mass of the first gauge KK state. 
The color coding of the curves corresponds to 5 = 5 n , , —5^ being represented by the 
top(bottom), middle (middle), bottom(top) curves on the right(left) handside of the figure. 
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achieved. We then reexamined the precision electroweak bounds on the lightest gauge KK 
mass or, equivalently, A„.. We found that for a reasonably wide range of parameters, the 
mass of the lightest gauge KK field can be as low as a few hundred GeV and hence will 
be visible at the LHC and possibly the Tevatron. The weakening of the constraints from 
precision data results in a significantly less severe bound on A n ; one can now naturally obtain 
A n < 10 TeV. This makes the RS framework with gauge fields in the bulk more favorable, in 
the context of a solution to the hierarchy problem, and facilitates the construction of more 
realistic models based on the RS proposal. In addition, for much of the 5 0jn parameter space, 
the gauge field and graviton KK excitations are nearly degenerate which may result in some 
interesting collider phenomenology. 

Given the observations in this paper, we expect that the addition of brane localized 
kinetic terms for other bulk fields, such as gravitons, could significantly alter the phenomenol- 
ogy of the RS model, and could also lead to the appearance of new features in the low energy 
theory. 
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Appendix 

Here, we will show that near special values of the BLKT coefficients 5q and 5 n a 
vanishing root develops in Eq. (jl5|) . signaling the appearance of an unphysical ghost state. 
We will also discuss the large \6 n \ limit of this equation. In the limit where x n —> 0, Eqs.(fT5j) 
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and (US} yield 

2 -4(^ + Sq + Ax c 7r) 

x « _( fcrcVr + 5 )( 2 ^ + l) + (l + 5 7r )- 1 J 

This equation is valid for x n <C 1, and for S n = — (S + kr c ii), it yields x 2 n = 0. Thus, for 
5-x < — (Sq + kr c 7r), an unphysical ghost state, a state with negative norm, appears in the 
spectrum, as implied by Eq. (J27|) . Note that in the original RS model, with Sq = S^ = 0, 
this ghost state does not appear, since kr c ir > 0. 

We now study Eq.(|3*Uj) in the limit 5 n — > ±oo, for the three cases S = 0, 5 n , —S n that 
were examined in the text. For So = 0, we have x 2 n — > 4/(2/cr c 7r — 1) as S^ — > ±oo. Since 
/cr c 7r Rj 35, the small x n approximation is valid and for S n — > +oo, the root is physical, 
as shown in Fig. la. Here, we see that the root for the lowest KK state asymptotes to a 
constant value of ~ 0.24 as S n becomes large and positive; when S^ is large and negative, 
the ghost root asymptotes to this same value. For S = S n , x 2 n — > 4/S^ as 5 n — > ±oo. This 
root is only physical for 5 n —> +oo, and approaches zero, as presented in Fig. |2K- In the case 
where S = —S^, we have x\ — > — (Akr c -n) / (2S 2 .) as S n — > ±oo, and there are no light physical 
modes in this limit, as seen in Fig. 2b. 

We would like to note that in the limit S n — > ±oo, with choices of Sq that yield 
\a n \ <C 1, Eq.(|I5|) is well-approximated by J\(x n ) = 0, for x n ^ 0. In this case, apart from 
the special light modes discussed above, the rest of the gauge field KK modes have masses 
that are approximately degenerate with those of the KK gravitons from the original RS model 
jS] . More generally, as can be seen from Figures and 121 a large range of parameters exist 
where this near degeneracy occurs. This will have important phenomenological implications 
at future colliders. 
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